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I. AMOS I ST PROGRAM 



A. HISTORY, PURPOSE, AND ORGANIZATION 

In the late 1960's the United States Army became aware that with the 
termination of the draft a source of physicians for the Armed Services 
would be lost. This loss would certainly have profound effects upon 
the heal tli services provided by the U.S. Army , especially in the areas 
of clinic care. In 1969 in order to help maintain adequate outpatient 
medical care in the face of manpower shortages, the Office of the 
Surgeon General formed the Automated Military Outpatient System (AMOS) 
Project. 

From the 'work of Project ADS lias come the ADS I ST Program. This 
program has emerged due to the ability of Project A v OS to record the 
logic used by physicians in their diagnosis and treatment of minor 
illnesses so that a nonprofessional medical person (ADSIST) with minimal 
training could safely treat certain types of medical problems. The 
ADSIST would then become a provider of health care, extending the 
capabilities of the attending physicians. The ADSIST Program is solely 
for ambulatory health care. 

The ATOSIST Program is made up of two parts called the Triage and 
the AIIC. In the Triage part of the program an ambulatory patient is 
screened by ai A'DSIST using a "Triage Note" (diagnostic decision tree 
questionnaire) in order to determine to which hospital clinic the patient 
should be sent. The Acute Minor Illness Clinic (AIIC) is staffed by 
ADSISTs and physicians (AMD SI ST MDs) who treat people screened for 
the AMIC by the Triage. The "Triage Note" will specify whether the 
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patient is to see an A'DSIST or an A'DSIST MD, but a patient may request 
to sec only an ASSIST MD if he desires. 

ADSISTs arc alv.'ays under the supervision of a physician, who sees 
patients sent directly from the Triage Section and patients referred 
from the various ADSISTs . An ADSIST may refer to a physician if he 
docs not feel sure of his diagnosis or is faced with an unfamiliar 
complaint. The physician may merely confer with the ADSIST verbally 
or lie may go to the A'DSIST's room to sec the patient himself. 

In order to make the program successful, the criteria for AMOSISTs 
are very stringent. Upon selection, these personnel are given two weeks 
didactic training at Fort Sam Houston, Texas, followed by about ten 
weeks on the job training as stated in Reference 1. Physicians attend 
the two weeks of didactic training and may then return to their facilities 
to help ready the clinic. 

At present the AMOS I ST Program lias been adopted in over 20 Army 
Medical Facilities. Adoption of this program at a particular facility 
is decided after studies have been done to ensure the program's feasibility. 

B. IMPLEMENTATION AT SILAS B. RAYS HOSPITAL 

In December, 1973, the A^DSIST Program was instituted at Silas B. 

Hays Hospital, Fort Ord, California. Figure 1 shows the physical layout 
of the AMIC. 

The Triage is physically separate from the AMIC and is located next 
to the patient records room. Patients who arc directed (Triaged) to the 
AMIC arc logged in at the AMIC reception desk and have necessary vital 
signs taken by a Nursing Assistant. The patient's folder is then placed 
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Amosist's Room 
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Reception 
Waiting Area 
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in either a physician or ADSIST rack depending on the "Triage Note" 
attached. 'Hie patient is seen when a server becomes available. Order 
of service is first in, first out. 

During a treatment, a chaperone may be necessary. Nursing Assistants 
or even AMOSISTs arc used for this. If referral to a physician is nec- 
essary, the AIDSIST will look for an available physician. He may have 
to wait if all physicians arc busy. When a physician is located, either 
a verbal referral will occur or the physician will go with the A DS 1ST 
to examine the patient. In either ease the effect is the same, namely 
that two persons arc occupied with one patient. 

A small percentage of patients come to the AMIC just for prescription 
refills. Their records are placed in a separate rack. Refill patients 
arc given a priority by the physicians and wait in a separate area from 
the other patients. Their service times arc normally very much shorter 
than those for other patients. 

Another small percentage of patients arc sent from the AMIC to the 
X-RAY Department or to the Laboratory. Upon their return to the AMIC, 
they are given priority by the .AMOS 1ST or physician who originally 
saw them. 

Figure 2 presents a block diagram of the system just described. It 
can be determined at this point that, with the exception of referrals 
from the AMOSISTs to physicians, the system can be viewed as two separate 
multi-channel queues. The analytical model will be discussed in another 
Section. 



11 



BLOCK DIAGRAM OF AMOSIST PROGRAM 
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II. AMIC DATA 



In order to develop an analytical model capable of analyzing and 
predicting the system behavior of the AMIC at Silas B.Hays Hospital it was 
necessary to establish a data base. To keep the modeling of the AMIC 
in a workable fora, all statistics related to prescription refills and 
patients sent from the AMIC to X-PAY or hAB were disregarded. Also, 
complete statistics were confined to weekdays when patient load was the 
heaviest. 



A. PERSONNEL STAFFING AND SCHEDULING 

During this study, there were four physicians available, one being 
responsible for the clinic's medical administration as well as seeing 
patients. There were 11 AUOSISTs available who either saw patients or 
assisted in the Triage. An additional non -medical enlisted person was 
the NCOIC. The support staff was made up of four Nursing Assistants 
and three receptionists. 

Tnc AdIC was open to receive patients from 0745 to 2330 hours, 
with active duty personnel having priority prior to 1000 hours . Tic day 
shift of personnel work until 1630 while the night shift worked from 
1530 to 2330 hours. During weekdays, three AMOS I ST 'IDs worked on the 
day shift and one at night . Since the AMIC was tasked with sending a 
physician to the stockade each morning, this meant that at most two 
physicians would be available during the mornings. 

On Monday , the heaviest day, seven to eight A'DSISTs would be 
scheduled to work the day shift and on the other weekdays, five to six. 

At night, two AMOSISTs were scheduled on all weekdays. 
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R. I’ATIENT ARRIVAL PATTERN 



In order to get an adequate estimate of the patient arrival pattern 
in the AMIC, log-in-sheets for February to May 1974 were reviewed. Each 
month's arrivals were recorded for each hour of the day and for each 
day of the week. These arrival figures were averaged together to give 
a standard week. Holidays and days affected by holidays were left out 
of the calculations since they were unusual in value and not representative 
of an average day. For periods that did not include a full hour period, 
the figures were adjusted to an hourly figure in order to make them 
comparative for statistical purposes. Table I summarizes the arrival 
figures for an average week. Each day's average was based on 17 weeks 
of data. 

Each day of the week exhibits approximately the same shape curve for 
arrival patterns over the length of the day and differs only in the 
average daily total number of arrivals. This allowed arrival figures to 
be pooled together giving one average weekday arrival pattern as depicted 
in Figure 3. Different patient volumes could then be achieved by merely 
moving the entire curve up or down to produce any desired patient volume. 

C. ADS I ST HD SERVICE TRIES 

There was no procedure within the AMIC to recover the length of time 
that an ADSIST HD spends with a patient or with an AMOSIST on a referral. 
To estimate these times, approximately 75 actual observations of patient 
encounters with. AMOSIST M's were made during randomly selected days over 
a three week period. The results of these observations showed that an 
A'DSIST MD spends approximately 13.55 minutes with his own patients and 
about 4.34 minutes of his time with an A'DSIST on a referral. 
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TABU: I 





ARRIVALS 


PER HOUR FOR 


AVERAGE 


LEEK 






MON. 


TIJES. 


LED. 


THUR. 


FRI. 


0745 - 0800* 


18.48 


15.24 


10.36 


9.88 


12.83 


0800 - 0900 


24.69 


23.50 


25.12 


22.00 


20.67 


0900 - 1000 


21.75 


20.81 


17.94 


16.65 


14.94 


1000 - 1100 


17.50 


13.56 


16.47 


13.76 


12.56 


1100 - 1200 


10.25 


9.00 


8.41 


S. 71 


9.67 


1200 - 1300 


7.88 


7.56 


6.12 


5.94 


6.72 


1300 - 1400 


11.31 


13.50 


12.47 


12.65 


11.06 


1400 - 1500 


10.19 


11.13 


9.71 


11.82 


9.44 


1500 - 1600 


9.06 


8.31 


7.65 


7.82 


8.50 


1600 - 1700 


8.06 


5.69 


6.24 


6.82 


6.11 


1700 - 1800 


7.31 


6.56 


5.12 


4.76 


4.33 


1800 - 1900 


5.44 


5.94 


6.71 


6.47 


4.39 


1900 - 2000 


6.25 


6.69 


5.71 


5.65 


4.67 


2000 - 2100 


5.13 


5.38 


4.47 


4.76 


4.33 


2100 - 2200 


3.56 


3.88 


3.06 


3.24 


2.83 


2200 - 2300 


2.00 


2.50 


2.55 


3.11 


3.22 


2300 - 2330'“ 


.50 


.76 


.88 


.88 


1.10 



* ADJUSTED TO HOURLY AVERAGE 



IS 



TYPICAL WEEKDAY ARRIVAL PATTERN 







tij ioi i >rjd stvaiiriv iwtw ;Jd\ti:lw 
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TIME. OP DAY 



D. ASSIST SERVICE TIMES 



Tlie service times for each A'DSIST were broken down into referral 
and non-referral categories just as for the physicians. The actual 
times spent with patients for each individual A'DSIST were taken from 
the Data Collection Sheets (computer generated questionnaire used by 
A'-DSIST to diagnose a patient's complaint) which are kept on file for 
each A'DSIST. These sheets also show by a physician's signature or an 
appropriate box being checked when a referral was made. After analysis 
of these actual patient encounter times, an average service time (time 
from start of one patient to start of next patient) was derived for 
each A'DSIST. In all but one case, at least 40 encounters were used to 
compute average "service times" for an A'DSIST. 

The results of this data indicated that the average service time for 
a non-refcrrcd patient was 17.96 minutes while the referred patient's 
service time was 23.85 minutes. The actual patient encounter times for 
the A'DSISTs compared quite favorably with results published in a recent 
ADS Project Progress Report (Reference 2). In addition, an average 
referral rate for an A'DSIST was found to be 30.93 o s. The above estimates 
were based on about 400 A'DSIST' s encounters. 

E. ‘ OTHER STATISTICS 

Since tlie AMIC is a system that generates patients to two different 
types of servers, it was necessary to find the percentage of people who 
were either triaged to see an A'DSIST MD or requested to sec an A'DSIST 
Ml). This type of data was not readily available from the records of 
the AMIC. For selected weeks the Data Collection Sheets for all A'DSISTs 
were totaled for certain days. This figure showed tlie total patients 
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triaged to an A'-DSIST who actually saw an A DS I ST on a given day. This 
was subtracted from the total patients entering the AMIC, as shown on 
the log-in-sheets, to give the number of patients seeing an ADSIST MD 
directly. This figure showed that approximately 55% of the patients 
were sent directly to the A'-DSIST MD versus 45% for the A'DSIST. This 
estimate originated from approximately 1000 cases. 

Since it was decided to start the analytical model at 0800 in the 
morning, an estimate of the average number of people in the waiting 
area prior to the first patient being seen was required. The arrival 
information previously gathered for February to May, 1974, was used to 
find an estimate of 11 people waiting to sec a server prior to the first 
sendee time. 

F. INTERPRETATION OF DATA 



It would be inaccurate to say at this point that the AMIC at Silas B. 
Hays Hospital can be fully represented with the aforementioned set of 
numbers. Any system that is composed of people and is dependent upon 
their interactions can probably never be completely and accurately 
modeled . 

The data gathered on the system were required by the analytical 
approach to the AMIC model. Also, the analytical model is a long run 
representation and average figures should not be expected to predict 
exactly what the situation will be like tomorrow. Some variation will 
almost always occur. 

Each system lias its own inherent characteristics that may not be 
readily discernible at first glance. The AMIC at Silas 15. Hays is no 
exception. The stockade commitment in the mornings requiring a physician 
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from the AMIC, the priority until 1000 for active duty personnel giving 
an uneven arrival pattern, and the makeup of the population that the AMIC 
serves, are things that may not be common to any other AMIC. Each 
system must be judged on its own performance in the face of its individual 
problems and not be compared using "raw data" to other systems which 
outwardly appear the same. 

The analytic model output to be presented pertain to Silas B. Hays 
liospital as it was operating during the time of this study and uses 
data taken at a particular period of time. No doubt as the system 
changes in time certain statistics may need revision. The greatest use 
of this study will be the chance to see the results of changes, either 
proposed or actual, and predict their effect upon the entire system. 

Table II summarizes the data that has been presented in this 
Section and which will be used as the data base for the analytic model. 
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TABLE II 



DATA SUMMARY 



ADSIST Ml): 

MEAN SERVICE (DIRECT) 13.55 min. 

MEAN SERVICE (REFERRAL) 4.34 min. 

ADS I ST: 

MEAN SERVICE (NON -REFERRAL) 17.96 min. 

MEAN SERVICE (REFERRAL) 23.83 min. 

REFERRAL RATE 30.93$ 



PROPORTION OF DIRECT .'ID PATIENTS 
MEAN ARRIVALS PRIOR TO FIRST SERVICE 
AVERAGE ARRIVALS PER WEEKDAY: 

MON. TOES. LED. 

156 149 145 



55.00$ 

11 



TOURS. 


FRI 


140 


127 



STAFFING: 

ATOSIST MD 4 
ATOSIST 11 
NURSING ASST . 4 
NCOIC 1 
RECEPTIONIST 3 
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III. Ull: AMIC AS A QUEUEING SYSTiH 



During the course of data gathering and research, the similarities 
of the Acute Minor Illness Clinic to a multi-channel queueing system 
became apparent. In the actual system the arrivals follow on inhomo- 
geneous pattern throughout the course of any day. Once in the system 
a patient is directed to an AMDS I ST or an AIDS I ST MD and then waits 
for the first available server. l\hen a server is free, the earliest 
arrival waiting is treated mid released. The major complication in 
modeling the actual system is the interference in the ADS I ST HD queue 
caused by referrals from ADSISTs. However , if the two queues (AMOSIST 
and ADSIST MD) can be isolated, two independent queueing problems can 
be formulated. 

The ADSIST referrals cause an interference in the AMOS I ST MD queue 
because the referred patient occupies an AMOS I ST MD with an additional 
service. This however, is the only effect because an AMOS I ST remains 
occupied while the referral takes place. Thus the problem can be 
simplified by determining an appropriate augmentation to the AMOSIST MD 
arrival stream to handle those patients referred by ADSISTs. These 
additional arrivals can be considered as "virtual" arrivals - i.e. 
persons who will eventually occupy the AMOSIST MD's time but who initially 
arrive to sec the .-ADSIST. 

The referral rate from ADSIST to ADSIST MD is taken as a percentage 
of the ADSIST arrival stream. The number of arrivals thus computed was 
then added to the ADSIST MD arrival stream. The original clinic system 
can then be viewed as taco distinct multi-channel queueing systems. 
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In order to use the data presented in Section II as an input to the 
analytic model, it was necessary to calculate certain new statistics. 

The hourly arrivals as well as a single average service time for an 
AM0S1ST and an ADSIST 'ID were computed. 

For a given hourly arrival figure, the product of this number and 
the average percentage seeing an A''DSIST gives the hourly arrivals for 
the AMOSIST. The hourly arrivals for an ADS I ST MD is the product of 
total hourly arrivals and the percentage seeing the physician directly 
added to the product of the referral rate and the hourly arrivals for 
an ADSIST. 

The analytic model recognizes only one average service time for an 
ADSIST or an A’TOSIST MD. For both servers, Section II gives two 
different types of service times. For the A'DSIST, it was referral 
and non-referral service times. The AMOSIST MD had either direct patient 
or referral service times. The percentage of each type of service given 
by an A'DSIST or A'DSIST MD was computed. Tie product of this percentage 
and the appropriate average service time added over types of service 
gives an effective average service time for the A'DSIST and ADSIST MD. 
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IV. Tin': MOPF.L 



A. GENERAL DESCRIPTION 

In order to analyze the two queueing systems, a mathematical model 
was developed. Using the iiihomogeneous arrival stream and server schedules 
as input, the model computes the probability distribution of system size 
throughout the day. The model then computes various measures of system 
performance at specified times of the day. The model treats the two 
queues as separate systems and during each day solves first the A'DSIST 
system and then the A’DSIST MD system. 

The significant problem in the formulation of the model was that 
all rates and schedules of servers were functions of time. There is no 
known analytical method of solution for the time dependent multi-server 
queueing system. Katadare and Kaufman (Reference 3) have advanced a 
Method of Successive Approximate Transient Solutions (MOSATS) for 
numerical analysis of a single channel time dependent system and much of 
their method has logical extensions to the multi -channel system. The 
approach calls for discretizing the time axis into segments over which 
the parameters remain constant. Then a transient solution is computed 
for each segment successively. The method was expanded to multiple 
channel queues and made efficient for real time computations for its 
application to this problem. The details are explained more fully below. 

B. MODEL ASSUMPTIONS 

The model makes several key simplifying assumptions . These will be 
listed and explained briefly here and referred to in subsequent, sections. 



(Al) Arrivals form an inhomogeneous Poisson Process. 

(A2) Services arc exponentially distributed and independent of 
the arrival stream. 

(A3) At most three arrivals and three departures can occur in any 
micro time step. A micro time step is a sub-interval of the periods 
in which parameters remain constant; it is used in the computation of 
state probabilities. 

(A4) No patient can both arrive and depart within the same micro 
time step. Patients leaving during a micro time step arc from those 
present at the beginning of the step. Tins assumption has the effect 
r of simplifying the calculations somewhat. 

C. J-DDhL Old: PAT I ON 



In order to adequately describe the operation of the model, a general 
description of the computations for solving a day's clinic operation is 
presented here. Mathematical amplifications arc presented with the 
appropriate steps where necessary. Output and interpretation are discussed 
later. 



1. Initial State Probability Vector 

The probability distribution of system size at the day's beginning 
is Poisson with a mean given by the user. Specifically, 

Pn(0) = ’ n = 

where P (t) is defined to be the probability that system size is n at 
time t raid p is the average number of early arrivals for the given day. 

To ensure that no probability mass was lost by limiting system size to 
N, the model places all remaining mass in the probability of largest 
system size. That is, 
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N-l 



P N (0) = 1.0 - 2 P n (0). 

H n=U 

In all applications this value was zero. System size (N) used in the 
model was 15 for computational efficiency. This restriction is not 
severe because in light of the two queue concept, it represents an over- 
all system size of 30. 

2. Time Step Computations 

Hie model next determines the first period of time during which 
all parameters (arrival rate and number of servers) arc constant. It 
should be noted that both arrival rates and server schedules are specified 
as step functions. Within this interval, subsequently referred to as the 
macro interval, time dependence is removed. Using (A3) and (A4) this 
macro interval is subdivided into micro intervals such that the assumptions 



arc valid. Specifically, 



1.25 



1 — ~ fJL • 

o - hr • 



1.25 

CM 



) 



where <5 is the length of the micro interval, A is the appropriate 
arrival rate, c is the number of servers and/* is the service rate per 
server. The length of the macro interval is then divided by 6 to 
determine the number of "steps" which will be necessary to compute the 
state probability vector at the end of the period. The constant (1.25) 
used in the 6 computation was experimentally determined as that figure 
which gave the best approximation to a correct steady state solution for 
a time independent system. 

5. Transition Probability Matrix 

i . ■ — — .... — . .. 

In order to determine the probability of entering a given state 
from an initial state in a £ step, a matrix of transition probabilities 
[Q] is next determined. . The entries in this matrix arc defined as 

Qi. n (t) = Pr (system size at t is n given system size at 
’ t-<£ was k) 
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Further define 



Q] ( • • (t) = Pr (i arrivals and j departures in 6 given 
* system size at t -6 was k). 

Then we have approximately that 

Q k>n (t) - £ 



1 ? J 

j >< 3 
k+i-j=n 



To determine the Q, two more definitions will be necessary: 

ic ,n 

-X(i) = Pr (i arrivals in a 5 step) i = 0,1, 2, 3 

Y(j,k) = Pr (j departures in a & step given k present 

at the beginning of the 6 step) j = 0,1, 2, 3. 

It is clear that the probability of an arrival is independent of system 

state and/or number of servers. In fact, by (Al) , the probability of 

i arrivals in 6 is 

i 

X (i) = (X4) c /i.' , i = 0,1, 2, 3 

By (A3) the number of departures in & is limited to three . This 
restriction defines three distinct cases for the probability of a 
departure based on system size and number of servers. These cases arc 
first discussed and then the combination of arrival probabilities and 
departure probabilities to determine the Qj, n (t) entries is discussed, 
a. System Size Less Than c 

An entering system size less than the number of servers 
implies that not all servers are occupied and a binomial probability of 
departure exists. Specifically, 

Y(j,k) =(^) (1.0-e-^ )j( c '^) k " j , j = 0,1, 2, 3 and j . 
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b. System Size Greater Than c+2 

When the beginning system size exceeds the number of servers 
by two or more, the system is saturated and no server becomes idle 
during the period for cases of at most three departures. The probability 
of a departure is therefore k- independent. Specifically, 

Y (j ,k) = (c/*6 )V C/<t * / j l , j = 0,1, 2, 3 and k>c+2. 

c. System Size Between c and c+2 

When the beginning system size is such that all servers are 
occupied but also night become idle at some point during the time 
interval , the exponential assumption (A2) must be interpreted properly 
to account for the possibility of idle servers. If, in fact, the servers 
become idle during & , the distribution of inter -departure time is 

2-Erlang followed by 3-Erlang, etc. depending on the specific number of 
servers becoming idle. Looking more closely at the specific cases 
involved, it can be shown that: 



for k = c+1, Y(0) = e 



-c/t^ 



Y(l) = (c/-«S) c 

2 , -(c-lU,. 

Y (2) = c (c v J ‘ ) (1-e -/ 4 de ) 



Y(3) = 1.0 - E Y ( j ) 

j=0 



and for k = c+2, 



Y(0) = e 



-CAxd 



Y(l) = (c /i6) e 



-C/ud 



Y(2) = (c/x<i) 2 e" C/tt<f /2! 



Y(3) = 1.0 - 






Y(j) 
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Finally, the entries in the transition probability matrix may be 

determined. By (A4) arrivals and departures are independent events and 

thus Q, . . = X(i) Y(j,k). It should be emphasized that the effects of 

K,i,j 

k (entering state) are in the Y(j) term. 

A necessary property of transition probability matrices in general 
is that the row elements sum to 1.0. That is, given any beginning state, 
the system must be in some state at the end of the interval. By (A3) 
the model discards any tail probability beyond three arrivals and/or 
three departures. To correct for this assumption and to ensure that 
correct row sums are obtained, the convention of placing all remaining 
probability mass in the X(3) and Y(3,k) terms was adopted, except in 
those cases where the entering state was less than three. In those cases 
the remaining probability was placed in the highest possible number of 
departures consistent with entering system size. (Recall (A4)). 

4. 'Hie Interval Computation 

Having determined the step transition probability matrix, it 
is necessary to show the computation used to determine follow-cn state 
probabilities. It can be readily seen that: 

for P(t) = (P Q (t), P^t),...,^)), 

P(t*< ) = P(t)K) 

r(t+ 2 <t ) = p(t+<J )[Q] = P(t)[Q][Q] = P(t)[0] 2 



P(t*i<J J » I'(tJtQ) 1 , 

where i is the previously computed number of 6 steps Contained in tiic 
macro interval under consideration. (See IV. C. 2 above). Vi thin the 
model, i is taken to be a power of two for computational efficiency. 
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5. Iterations 



The results of the preceding computations are state probability 
vectors for the end of each macro interval. This state probability 
vector replaces the initial probability vector (1) and the computations 
(2 through 4) are repeated iteratively until a predetermined end of day 
is reached. The model then returns to the second arrival rate and server 
schedule (/CDS I ST MD) and repeats the entire procedure until the pre- 
determined end of day is again reached. This then represents the end of 
one day's clinic operations. 

D. STATISTICS 

1. Identified 

It was initially recognized that probabilities of system size 
and delay would be necessary if the model were to be of any practical 
value to the user. After considering the myriad of available interpretive 
tools, the list was ultimately narrowed to eight distinct measures of 
the system which are determined at each time the user desires a report. 
These statistics are defined and discussed separately in the following 
sections. 



a. Expected System Size is computed as 



£ nP (t). 
n=0 



' b. Standard Deviation of System Size is computed as 




c. The probability a server is busy at t is 



5 = (j; nP (t)+c P (t))/c. 

n=n 1 n =<’ n 



n=0 n=c 
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It should be noted that this probability is analagous to the standard 
utilization determination of queueing theory and is hereafter referred 
to as utilization (ori>). 

d. Tire expected size of the queue at t is 




nP (t) 
i-r J 



c-1 



nl> (t) 



n 



N 

c ZT P n (t) 

n=c n 



e. The standard deviation of queue size at t is 

Vt (n-c) 2 P n (t). 
n=c 

f. The probability an arrival at t encounters no delay is 

^ P ft') 
n=TT n ltJ * 

g. The probability an arrival at t encounters delay in excess 
of d minutes is 

N 



£ \ I" (gl c ^W<’J K | r (t). 

n=c \ k=0 K - 



h. The expected delay for an arrival at t is 




0 



2. Interpretation 

It is worthwhile to note that because of the augmentation of the 
ADSIST 1 ID arrival rate it is necessary'' to adjust certain of tire above 
statistics to remove the effect of the "virtual" arrivals. Statistics 
a,b,d and c are therefore interpreted as follows: a and d: the 
expectation of system and queue size are adjusted by a factor of p where 
p represents the percentage of the arrival stream which are actual 
ADSIST MD arrivals; b and e: the standard deviations are computed as 






fpqh(X) + p Y(X) 
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where p is as above, q is (1-p) mid X represents the computed and 
unadjusted figures for b and e. 

E. IBM 360/67 APPLICATION 

1. General 

The above model has been programmed in G-Lcvcl Standard Fortran 
for implementation on the IBM 360/67 under OS and for the CP/CMS time 
sharing system at Naval Postgraduate School. The general flow chart of 
the model is included as Figure 4. The model successfully reflects the 
actual system operation of the AMIC at Silas B. Hays Hospital. It lias 
been incorporated in an interactive time sharing mode and through the 
use of portable terminal equipment has been denonsti'atcd to hospital and 
clinic administrators. 

2. Input 

The model is sufficiently general to allow user control of the 
following input parameters: 

a. Volume. 

The user may specify any number of arrivals desired in the 
clinic for a given day. The model interprets this volume as a scale 
factor for the standard arrival rate step function. 

b. Service Time 

The user may vary the mean service time for either type of 
server and also for referral and non-referral service. 

c. Schedules 

'Hie user may specify schedules for each type of server and 
may have the number of each type server on duty change as many as 16 
times throughout the day. 



31 



FIGURE 4 



FLOWCHART 
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d. Percent of Arrival Stream 



'Hie user may specify that percentage of the arrival stream 
which sees the doctor directly. 

e. Referral Rate 

The user may specify that percentage of patients which will 
ultimately see a doctor through referral from an A N OS I ST . 

f. Early Arrivals 

11 ic user may specify the number of arrivals before the system 

opens. 

g. Reports Desired 

The user may specify up to 16 times throughout the day at 
which he would like statistics computed and displayed. 

h. Delay Cutoff 

'Die user may specify the number of minutes for which he 
would like the probability of delay in excess of that figure computed. 

(the "d" in IV.D.l.g above). 

3. Output 

Although eight statistics arc generated within the model, it is 
apparent that not all of these have meaning for the non statistically 
oriented user. For this reason only five of the eight arc normally 
displayed in interactive user sessions. These five arc tire of the 
report, utilization, expected size of the queue, expected delay encountered 
and the probability of delay in excess of the user supplied delay cutoff. 
The entire list of statistics is printed off-line. 

Tables 3 and 4 show the statistics generated for a typical 
weekday for A'DSIST and A'DSIST MP respectively. The server schedules 
arc also shown. The normal method of gathering statistics is to compute 



33 



at one minute past each hour to capture conditions which will be in 
force for the coming period. That is, because most changes in schedule 
and arrival rate step functions occur on the hour, the statistics are 
taken during the following period so that the conditions which will 
attempt to clear the system will be in force. 
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TABLE III 

AMOS I ST QUEUE STATISTICS FOR A TYPICAL WEEKDAY 
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AMOS I ST MD QUEUE STATISTICS FOR A TYPICAL WEEKDAY 
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V. CONCLUSIONS 



A. VALIDATION AND SENSITIVITY 

The computer model was validated by forcing "steady state" solutions 
to be achieved. Tills can be easily done by assuming a constant moan 
arrival rate, a constant number of servers and a constant mean service 
rate. In so doing, it became apparent that the model is sensitive to 
several parameters. 

1. Early Arrivals 

The effect of arrivals before the system is open for business 
cannot be overemphasized. The effect of the early arrivals is intuitive 
in that if the system begins congested, then it must spend a good deal 
of time alleviating this congestion. Arrivals during the initial periods 
will normally experience much longer delays than if no patients arrived 
early. Within the model, early arrivals is an input parameter. 

2. ASSIST MD Schedule 

As the output tables indicate, during normal operation the major 
system congestion occurs in the A^DSIST MD queue. This is attributable 
to the percentage of the arrival stream which is directed to the AMOS 1ST MD 
and to the augmentation of this arrival stream by referrals from the 
AIDS I ST as well as the number of .ADSIST MDs available. This model is 
very sensitive to the ANOSIST MD schedule and a large expected delay can 
be seen when few servers arc scheduled during peak periods. In actual 
clinic operation the A'DSIST MD arrival rate is affected by those patients 
who request to see an A'DSIST MD and also by those ADSISTs who refer a 
higher percentage of their patients. The percentage of arrivals going 
directly to the A'DSIST MD and the referral rate from AMOSISTs to A'DSIST 
MDs arc inputs to the model. 
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B. StMIARY 



' Ihc Acute Minor Illness Clinic at Silas B. Hays Hospital can be 
successfully modelled os a multi -channel, time dependent queueing system 
and approximate solutions to desired queue and delay statistics can be 
found using numerical methods. 

The computerized version of the model is capable of handling much 
larger systems and could be expanded to handle more types of servers. 

As previously mentioned, the present version with system size limited 
to 15 is used interactively with the user supplying input and seeing 
output statistics within 30 seconds. 

A general method of solution for multi-channel, time dependent 
queueing systems has been presented. This model could find ready appli- 
cation to any number of similar queueing systems presently unsolved. 

The procedure presented is sufficiently general so that a wide variety 
of, and variance in, input parameters may be easily handled thereby 
enhancing the model's applicability. 
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